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ABSTRACT 



We report the discovery made using the Austraha Telescope Compact Array 
of a remarkable string of radio emission towards IRAS 16562-3959, a luminous 
infrared source with a bolometric luminosity of 7.0 x 10^ L© . The radio emission 
arises from a compact, bright central component, two inner lobes, which are 
separated by about 7" and symmetrically offset from the central source, and two 
outer lobes which are separated by about 45". The emission from the central 
object has a spectral index between 1.4 and 8.6 GHz of 0.85 ± 0.15, consistent 
with free- free emission from a thermal jet. The radio emission from the lobes 
have spectral indices in the range characteristic of thermal emission. We suggest 
that the emission from the lobes arises in shocks resulting from the interaction 
of a coUimated wind with the surrounding medium. The radio string is located 
within a massive dense molecular core, and is associated with extended green 
emission (Spitzer 3-color), Herbig-Haro type emission (2MASS K^-band) and 
OH maser sites - all phenomena readily observed towards sites of massive star 
formation. We conclude that the massive core hosts a high-mass star in an early 
stage of evolution in which it is undergoing the ejection of a powerful coUimated 
stellar wind, showing that jets found in the formation of low-mass stars are also 
produced in high-mass stars. 

Subject headings: ISM: (IRAS 16562-3959) — ISM: jets and outflows — radio 
continuum: stars — stars: formation 
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INTRODUCTION 



The determination whether massive stars are formed via accretion or via merging 
processes is one of the main observational challenges in the field of star formation. If massive 
0-type stars are formed by an accretion process similar to that inferred for low-mass stars, 
then we expect that circumstellar disks and jets will be present in their earliest stages of 
evolution. 

To date there are only a handful of massive young stellar objects (YSO) known to 
be associated with highly collimated jets and/or Herbig-Haro (HH) objects. All except 
one have luminosities smaller than 2 x 10^ Lq correspon ding to that of a BO ZAMS 



star. They include IRAS 18162-2048 (L ~ 1.7 x 10^ L 



HW2 fL ~ 1 X 10^ ; 



Cesaroni et al. 



1999 



Rodriguez et al. 



1 



Marti et al. 



19931 ): Cepheus A 



19941 ): IRAS 20126+4104 (L ~ 13 x 10^ Lg, : 



1997h: G192 . 16-3.82 (L ~ 3 x 10=^ Lq ; 



Shepherd et al. 



Shepherd et al. 



1998 



Devine et al. 



20011 ): and W75N, which cont ains several molecular and HH outflows 



powered by at least four late to early-B protostars (IShepherd et al. 



2OO3I ). There is only 



Rodriguez et al. 



Beuther fc Walsh 



Garav et al. 


2003; 


Rodriguez et al. 


2005; 


Brooks et al. 


2007 



2008h. Two other luminous YSOs, IRAS 180 89-1732 (L ~ 3.2 x 10 ^ L© ; 



Bronfman et al. 



20081 ). are 



20081 ) and G331.51-0.10 (L ~ 1 x 10^ Lq ; 
associated with radio continuum sources with spectral indices characteristic of collimated 
stellar winds, but the angular resolution of the observations is insufficient to resolve the 
jet/fiow morphology. It is not clear whether the lack of young massive stars with spectral 
types earlier than BO ZAMS associated with jets and/or disks is an intrinsic property of 
the most massive stars or due to observational disadvantages - massive stars are rarer and 
their evolutionary time scales are much shorter than those of low mass stars. 



Currently we are carrying out a systematic search for jets towards massive YSOs to 
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assess whether or not they are a common phenomena (Guzman 2011, in preparation). In 
this paper we report the discovery, made using the Austraha Telescope Compact Array 
(ATCA), of a string of radio continuum emission associated with IRAS 16562—3959. 
Assuming that this so urce is located at the distance of 1.6 kpc (Vlsr = —12.6 km s~^. 



Urquhart et al 



, the total 



IRAS fluxes (see iCasoli et al. 



ar-infrared luminosity, computed using the observed 
1986[ ). is ~ 7.0 X lO'' Lq . The target was selected from a 
list of sources compiled by us with IRAS luminosities in excess of 2 x 10^ Lq and radio 
emission much weaker than that expected from the total FIR luminosity. The expectation 
is that objects with these characteristics are high-mass objects in the pre-UC H II region 



phase, in w^ 
(see Fig. 6, 



lich the weak radi o emission is most likely arising from stellar wind phenomena 



Hoareeta 



(G345.4938+ 01.4677, 



2007). IRAS 



Urquhart et al. 



with ATCA (Urquhart et al 



6562— 3959 was observed as part of the RMS Survey 



2008al ). The emission detected at 4.8 and 8.6 GHz 



20071 ) was weaker than expected. Moreover, the positive radio 



spectral index between these two frequencies and the hint of a string of multiple sources 
gave us strong indications that this source could be an ionized jet associated with a massive 
young stellar object (MYSO). 

In this paper we present new ATCA data taken at 4 frequencies (1.4, 2.4, 4.8 and 
8.6 GHz). Based on these data, we propose that the string of five radio sources detected 
towards IRAS 16562—3959 is composed of a central thermal jet, plus inner and outer lobes 
of shock ionized gas resulting from the interaction of the highly coUimated stellar wind with 
the surrounding medium. 



^The two-fold distance ambiguity was resolved in 



Faundez et al. 



(120041 ) 
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2. OBSERVATIONS 

The radio continuum observations were made using the Austraha Telescope Compact 
Array (ATCA) I during 2008 June, October and 2009 February. We used the 1.5B, 1.5C 
and 6. OA configurations, utihzing all six antennas and covering east-west baselines from 
30 m to 5.9 km. Observations were made at four frequencies: 1.384, 2.368, 4.800 and 8.640 
GHz, each with a bandwidth of 128 MHz, full Stokes. Throughout this work we will refer 
to these frequencies as 1.4, 2.4, 4.8 and 8.6 GHz, respectively. The phase center of the array 
was RA= 16^59"41.61^ DEC = -40°03'43.4" (J2000). The total integration time at each 
frequency was about 180 minutes, obtained from 10-minute scans taken over a wide range of 
hour angles to provide good (u,v) coverage. The calibrator PKS 1740-517 was observed for 
3 min. before and after every on-source scan in order to correct the amplitude and phase of 
the interferometer data for atmospheric and instrumental effects as well as to calibrate the 
bandpass. 

The flux density was calibrated by observing PKS 1934-638 (3C84) for which values 



of 14.95, 11.59, 5.83, and 2.84 Jy were adopted at 1.4, 2.4, 4.8, and 8.6 GHz, respect 



Standard calibration and data reduction were performed using MIRIAD (iSault et al. 



ively. 



19951). 



Maps were made by Fourier transformation of the uniformly weighted interferometer data. 
The noise level achieved in the images are 0.38, 0.21, 0.096, and 0.070 mJy beam~^ and the 
synthesized (FWHM) beams obtained were 10.23" x 5.71", 5.97" x 3.27", 2.81" x 1.83", and 
1.62" X 1.01", at the frequencies of 1.4, 2.4, 4.8, and 8.6 GHz, respectively. We estimate the 
flux uncertainty to be approximately 20%. 



^The Australia Telescope Compact Array is funded by the Commonwealth of Australia 
for operation as a National Facility managed by CSIRO 
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3. RESULTS 

Figure [1] shows maps of the radio continuum emission from IRAS 16562—3959 at the 
four frequencies observed with ATCA. The maps at the lower frequencies (1.4 and 2.4 GHz) 
show three sources roughly lying in a linear structure along a direction with P.A. of 99.1°. 
The outermost sources which we will call Outer-East and Outer- West, labeled 0-E and 
0-W in Fig. [H respectively, are separated by an angular distance of ~ 45" (corresponding 
to 0.35 pc at 1.6 kpc). The maps at the higher frequencies (4.8 and 8.6 GHz) show that 
the central object seen at low frequencies is resolved into three components aligned along 
a direction with P.A. 98.9°. The external components of this triplet, which we will call 
Inner-East and Inner- West components (labeled I-E and I-W) are symmetrically located 
in opposite directions from the bright central source (labeled C), and separated by ~ 7". 
Throughout this work, we will use the abbreviations C, I-E, I-W, 0-E, and 0-W for the 
Central, Inner-East, Inner- West, Outer-East, and Outer- West components, respectively. 

We note that the five radio components are not exactly aligned, exhibiting a small 
bending. The lines joining the east components and the central source and the west 
components and the central source form an angle of 168°. We also note that in the 8.6 GHz 
image there are three small knots located about 4" north from the Outer-East lobe, with 
peak fluxes of 0.6 mJy. The knots are well aligned with the Inner and Central components 
suggesting that they could be part of a string of emission, but more sensitive observations 
are needed to confirm this or disprove as artifacts of the data reduction. We will not discuss 
these components further in this work. 

The position and flux densities of all radio components detected towards IRAS 
16562—3959 are given in Table [1] For all sources the flux densities were determined 
from Gaussian fittings using the MIRIAD task IMFIT, with the exception of the outer 
components at 4.8 and 8.6 GHz. For these, the flux density was measured over the 
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3(7 contour level. To estimate the flux densities of the C, Inner-East and Inner- West 
components at 1.4 and 2.4 GHz we fitted the observed morphology with three unresolved 
Gaussians centered at the positions determined from the 8.6 GHz observations. 

Figure [2] shows the SED of the integrated radio continuum emission from the five 
cor nponents. Also plotted in the spectra of the central component are the values reported 



by 



Urquhart et al. 



(120071 ) at 4.8 and 8.6 GHz of 4.8 and 12.5 mJy, respectively. We find 
that the radio continuum spectra of component C is well fitted by a power law spectrum 
{Su oc z/", where Si, is the flux density at the frequency z/), with an spectral index a of 
0.85 ± 0.15. For the other components the data were fitted using a model of thermal 
emission from a homogeneous region of ionized gas. The parameters EM and Q, emission 
measure and solid angle, are given in each panel. 

4. DISCUSSION 



4.1. The nature of the radio sources 



The spectral index of the radio continuum emission from component C, of 0.85 ± 0.15, 
indicates free-free emission arising from a thermal jet. Theoretical calculations show that 
collimated stellar winds can have spectral indices in the ra nge 0.25 to 1.1 depending on 



19861). They also 



the radial dependence of the physical quantities of the jet ( iReynoldsl 
predict that the angular size depends with frequency. However, the angular resolution of 
our observations does not allow us to investigate this dependence and observations with 
higher angular resolution are needed to resolve the jet. 

Assuming that the Central source is a bipolar, pr e ssure confined jet, which has an 



spectral index of 0.84, Eqns.(13) and (19) of iReynoldd ( 1l986l ) yields, for a distance of 1.6 



kpc and an observed flux density at 8.6 GHz of 12.1 mJy, the following constraints on the 
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jet physical parameters: 



106 



cm 



1.32 



2.8 



10-3pc 



-0.9 



zy,„ ^^ 0.818 ^ ^ sin z ^ 

sin 45' 



lOGHz/ V0.2 



r 



103 km s-V VlOGHz/ V 0.2 



sin 45° 



-1/4 



(2) 



10-6 M© yr-i 

where njet(?") is the number density of the jet at a distance r from the MYSO, z/^ is the 
turnover frequency, 6 the opening angle at the base of the jet, i the inclination angle with 
respect to the line of sight, the wind velocity and My^ is the mass loss rate. Most of these 
parameters are unknown. In deriving these expressions, in addition to the observational 
constraint, we assumed an ionization fraction of 1 and a temperature of 8000 K at the base 
of the jet; and a mean particle mass of 1.3 mn per hydrogen atom. 

The degree of collimation of the jet source can be estimated from the size of the lobes 
and their distance to the central component. Using the data obtained for the East and West 
outer lobes we derive that 6 ~ 0.2 radians . Assuming th a t Vm is 500 km s" ^ , a value typi cal 



of iets associated with lumi nous objects (lAnglada 



2006 



Rodriguez et al. 



1996 



Marti et al. 



1998 



Curiel et al 



20081 ). then the constraint equation ([2]) implies that the mass loss 



rate is ~ 1.4 x lO'^ Mq yr'^ 

If the total luminosity of IRAS 16562—3959 is produced by a ZAMS star it would 
corre spond to an 08 star, which emits a rate of ionizing photons of 2.2 x 10^^ s~^ ( jPanagia 



19731 ). Embedded in a constant density medium, this star would generate an H II region 
with a flux density of ~ 8.4 Jy at optically thin radio frequencies, far in excess of the 
observed value of ~ 10 mJy. We suggest that the weak radio emission from the central 
source is a consequence of IRAS 16562—3959 undergoing an intense accretion phase, with 
the central object still being in the pre-UCH II region sequence phase. The high-mass 
accretion rate of the infalling material forbids the development of a sizeable H II region 



( lYorke 



1979 



Walmsley 



19951 ). and the free- free emission from the ionized material is 



- 9 - 



considerably lowered at centimeter wavelengths. 

The spectra of the radio emission from the lobes are rather fiat, indicative of optically 
thin free-free emission. We suggest that the radio emission from the lobes arises in shocks 
resulting from the interaction of a coUimated stellar wind with the surrounding medium. 
The dotted lines in Fig. [2] correspond to fits of the observed spectra using a model of 
thermal emission from a homogeneous region of ionized gas. From these fits we derive 
emission measures of ~ 1 x 10^ pc cm^^ towards the inner lobes and ~ 1 x 10^ pc cm~^ 
towards the outer lobes. For all lobes the emission is optically thin within the range of 
observed frequencies. The solid angles were derived from the radio images giving 3.7, 
0.8, 25 and 15 arcsec^ for the Inner-East, Inner- West, Outer-East and Outer- West lobes, 
respectively. The derived electron densities within the lobes are about ~ 10^ cm~^ for the 
outer lobes and ~ 10^ cm~^ for the inner lobes. 



4.2. Characteristics of the exciting source and its environment 



The radio string is associated with a mass ive 1.0 x 10 ^ M^;^] and dense (1.2 x 10 



cm ^) core detected at 250 GHz (1.2 mm) by 



Faundez et al 



(120041 ) ■ with the Central 



component offset from the millimeter continuum peak by ~13", equivalent to 0.1 pc at 1.6 
kpc. 

Figure [3] shows that associated with the dust core there is strong diffuse band 
emission, centered near the Central source and aligned along the jet axis. This emission 
extends ~ 1.5' in the East and West directions. There is also diffuse Kg band emission 
arising from two regions located ~ 3' towards the North and South of the Central 
source. They are roughly symmetrically displaced from the Central source but aligned 
in a perpendicular direction to the jet axis. The strong band emission is likely to 
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be produced by excited H2 2.12 fim emission arising from shocked gas. Moreover, the 
arc-hke morphology of the northern band emission feature located approximately at 
RA=16''59'"42^ DEC= -40°00'13" (J2000) is also consistent with HH phenomena. We 
note that the North-South structure is also seen as diffuse mid-infrared emission in the four 
MSX bands (8.28, 12.13, 14.65, and 21.34 /xm ), indicating also a strong contribution from 
warm dust. 

Figure H] shows a zoom in of the 8 fim emission from Spitzer and band emission 
from 2MASS towards the central part of the core. There is an extended "V" shape feature 
extending along the E ast jet axis seen in both band and 8 fim . In addition, there 



are three OH masers ( Caswell 



1998 



2OO4J ): one is associated with the central source, a 
second is associated with a 8 /im source, unseen in 2MASS, located ~ 15" northeast of 
the central source and the third is located close to the East jet axis at about 2" from 
the Outer-East lobe. Th eir velocities are close to the radial velocity of the ambient cloud 



(V 



LSR 



-12.6 km s-\ 



Urquhart et al. 



2008bl ). except for the eastern OH maser, which 



has a radial velocity of —24.5 km s ^. Infrared continuum emission at 10.4 am is shown in 



Mottram et al. 



(120071 ). It shows 



Fig. [5l which displays an imagcl obtained with TIMMI2 by 
intense emission associated with the Central source and more diffuse emission associated 
with the Inner-East radio lobe. There is also an extended 4.5 fim emission enhancement 
associated to the central source (see Fig. [6]). These "green fuzzies", as seen green in a 
three color IRAC images, have been related to shocked gas in protostellar environments 



(jChambers et al. 



20091 ). All of the phenomena described above is readily associated with 
high-mass star formation and therefore support the notion that IRAS 16562-3959 is a 
high-mass young stellar object with energetic outflow activity. 



While the East- West K,, band structure can be related directly to the radio jet, the 



■^Downloaded from 



http:/ /www.ast. feeds. ac.uk/RMS/ 
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identification of tlie energy source of tfie Nortfi-Soutli band feature is less clear. One 
possibility is that the North-South feature is driven by the Central jet source, which will 
imply a re-orientation of the jet in almost 90°. Models of close stellar encounters do n ot, 
however, predict such large precession of the outflow axis (e.g. 



Moeckel &: Ballv 



20061 ) 



and thus this possibility appears unlikely. Since massive YSOs are known to be formed 
in clusters, the most plausible explanation is that a yet unidentified nearby source is the 
exciting source of the North-South K^-band emission. A possible candidate is the YSO, 
located 16" north-east of the Central component, identified from the Spitzer data. As 
pointed out previously, this object is associated to an OH maser with a radial velocity 
similar to that of the ambient cloud. This and the Central source belong thus to the same 



star forming region, as suggested by 



Caswell! (|l998h . 



Figure [7] shows the spectral energy distribution (SEP) of IRAS 16562-3959 from 3.6 



fim to 1200 fim , including flux densities at 1.2 millimeter ( iFaundez et al. 



2004h : at 12, 25, 



60, and 100 fim obtained from the IRAS database; at 8.3, 12.1, 14.7, and 21.3 fim obtained 
from the Midcour se Space Experiment (MSX) Survey of the Galactic Plane database 



(^Price et al. 



2001 



MSX source G345. 4938+01. 4677), and at 3.6 and 4.5 /im obtained from 
Spitzer-IRAC images using an aperture of ~ 1'. The IRAC images at 5.8 and 8.0 /im data 
were saturated and an adequate flux estimation was not possible. In the frequency range 
covered by the SED the emission is mainly due to thermal dust emission. The spectral 
energy distribution was analy zed by fitting model SEDs using a large grid of precomputed 



models (IRobitaille et al. 



120071 ). The continuous line in Fig. [7] presents the result of the best 
fit. The derived parameters from the fit are given in Tab. |21 The SED fitting indicates that 
the energy source corresponds to a deeply embedded YSO, and rules out the possibility of 
being a heavily extincted star. 



The mass of the host envelope derived from the fitting is ~ 1700 Mq. However, 
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the model of iRobitaille et al.l ( 120071 ) does not consi der the increase of absorptio n at mm 



19941 ) and in 



frequencies due to ice coatings in the dust grains ( lOssenkopf fc Henningl 
consequence is probably overestimating the mass of the envelope by a factor 2^3. 
Therefore, the mass of the core derived from the 1.2 mm observations of 910 Mq is in fair 



Faundez et al. 



(I2OOJ) reported a mass of 



agreement with the SED fitting. We note that 
1500 Mq, but a new analysis of the data shows that for this particular source an erroneous 
calibration factor was applied, which produced an overestimation of the fiux by a factor of 
1.76 (Lopez et al. 2010 in preparation). 

The fitting also indicates that the central YSO has a mass of ~ 15 Mq and a luminosity 
of 5.8 X 10^ Lq , that the envelope is undergoing an intense accretion phase with a rate 
of 3.4 X 10^'^ Mq yr~^ and that the accretion rate onto the central object is 5.5 x 10~^ 



Mq yr ^. The derived ste 



star (ISternberg et al. 



lar parameters are simila r to those of a B0~09.5 main sequence 



20031) and a 08.5 ZAMS star flPanaria 



I973I ). We also note that the 



derived infall accretion rate is greater than ~ 1.0 x 10 ^ Mrr. yr ^ , the infall rate needed to 



quench the development of a sizeable H II region (IWalmsley 



19951). 



The linear morphology and spectral characteristics of the radio sources suggest that 
the string is physically associated with a highly collimated wind arising from a young, 
massive YSO with the lobes tracing the interaction between the jet and the ambient cloud. 
We postulate that the lobes radio emission corresponds to thermal emission coming from 
shocked ionized gas situated in the working surfaces of the Central jet source. Particularly, 
we rule out the possibility that the radio emission has an extragalactic origin - the 
probability of finding an extragalactic source at 5 GHz with a fiux density above 5 mJy 



within 20" from the peak of the 1.2 mm dust emission is 3.5 x 10 ^ ( iFomalont et al. 



199l|). 



- 13 - 



4.3. Jet and Shock parameters 



In this section we make use of the observed quantities of the lobe emission to derive 
parameters of the jet-ambient medium interaction. Considering that the free-free emission 
from the lobes arises from a shock wave formed at the head of the jet, th en the emission 



measure, EM, is given in terms of shock parameters by ( ICuriel et al. 

1.68 



19931) 



EM 



1.39 



0.8 



(3) 



106 pc cm-6y VlQS cm-V \^100 km s^V VlO'^K 

where n^, Vg and Te are the particle ambient density, the shock velocity and the temperature 
of the shock-ionized gas. If the ambient density is known, this expression can be used to 
estimate the shock velocity. 

To estima te the ambient density at the position of the lobes we re-analyzed the 1.2-mm 
data taken by iFaundez et al.l toward IRAS 16562—3959 and fitted the observed radial 
intensity profiles from the core region using a mod el in which the density and temperature 



follow power-law radial distributions (lAdams 



199ll ). The best fit to the image was obtained 
with the following molecular gas density and temperature profiles: na{r) = nQ{0.1 pc/r)^-^ 
and T(r) = To(0.1 pc/r)°'^, where uq = 1.4 x 10^ cm~^ and Tg = 50 K. We assumed a 
gas to dust mass ratio of 100 and a mean molecular mass of 2.3 mn- The derived density 
profile implies that the mass of the core, within a radius of ~ 90", is 790 Mq, and that the 
ambient density at the position of the East and West inner lobes - displaced from the core 
center by 10" and 18"- are 2.3 x 10^ and 7.6 x 10^ cm~^, respectively. Using these densities 
and the emission measures computed in §3l we derive from Eq. (|3]) that the shock velocities 
at the East and West inner lobes are 119 and 91 km s~^, respectively. The same analysis 
applied to the East and West outer lobes gives shock velocities of ~24 and ~66 km s~^, 
respectively, consistent with the expectation that shock velocity should decrease away from 
the jet source. 
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Assuming a simple one dimensional isothermal shock model (e.g. iMasson fc Chernin 



19931 ) and that the momentum flux imparted to the shock ionized material by the jet equals 



the ram pressure of the ambient material, then 

{V,-V,fp, = V^Pa, (4) 

where pj and Vj are the mass density and velocity of the jet, respectively, and pa is the 
ambient gas mass density. If the shock velocity, ambient density and jet density are known 
this expression can be used to estimate the jet velocity. The density of the jet at the 
positions of the inner lobes can be estimated from Eq. ([1]). The lobes are displaced from 
the jet by 3", which translates into a physical distance of 0.033 pc, assuming an inclination 
of 45°. Using 9={].2 and assuming 2/^= 10 GHz, then the density of the jet at 0.033 pc is 
1.2 X 10~^^ gr cm~^, about ~ 4 times lower than the ambient density. Using the values 
derived for the East and West inner lobes of the shock velocity, ambient density and jet 
density we determine from Eq. (jlj) jet velocities of 360 and 280 km s~^, respectively. 

We further estimate that the most recent ejection, giving rise to the inner lobes, took 
place between 90 and 120 years ago. Moreover, from Eq. ([2]) and taking again 10 GHz as 
the turnover frequency, we derive mass loss rates of 1.0 and 0.79 xlO~^M0 yr~^, which in 
combination with the ejecta velocities give momentum rates of 3.7 x 10~^ and 2.2 x 10""^ 
Mq km s^^ yr^^ for the East and West halves of the bipolar jet, respectively. 



5. SUMMARY 

We made radio continuum observations at 1.4, 2.4, 4.8 and 8.6 GHz, using ATCA, 
toward IRAS 16562-3959, a luminous object (L ~ 7 x 10"^ Lq ) thought to be a massive 
star forming region in an early stage of evolution. The main results and conclusions are 
summarized as follows: 
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1. The radio continuum observations show the presence of a remarkable string of 
radio emission, consisting of a compact, bright central component, two inner lobes, 
separated by about 7" and symmetrically located from the central source, and two 
outer lobes, separated by about 45". 

2. The emission from the central object has a spectral index between 1.4 and 8.6 GHz of 
0.85 ± 0.15, consistent with free-free emission from a thermal jet. Assuming that the 
jet corresponds to a bipolar pressure confined wind with an aperture angle 9 — 0.2, 
we estimate that the jet has a total mass loss rate of ~ 2 x 10~^ yr~^. 

3. The radio emission from the lobes have spectral indices of typically —0.1, characteristic 
of thermal emission. We suggest that the emission from the lobes arises in shocks 
resulting from the interaction of a coUimated wind with the surrounding medium. 

4. The string is found projected towards a massive and dense core and associated with 
several indicators of massive star formation and shock gas tracers. The jet is located 
near the peak position of the dust emission. We conclude that the massive core hosts 
a high-mass star in an early stage of evolution in which it is undergoing the ejection 
of a powerful coUimated stellar wind, showing that jets found in the formation of 
low-mass stars are also produced in high-mass stars. 

The authors gratefully acknowledge support from CONIC YT through projects 
FONDAP No. 15010003 and BASAL PFB-06. This pubhcation made use of the 
GLIMPSE-Spitzer database. This paper also made use of information from the Red MSX 
Source survey database at www.ast.leeds.ac.uk/RMS which was constructed with support 
from the Science and Technology Facilities Council of the UK. 
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Table 1. OBSERVED PARAMETERS OF RADIO SOURCES 



Source 8.6 GHz peak position Flux density (mJy) 

a(2000) (5(2000) 1.4 GHz 2.4 GHz 4.8 GHz 8.6 GHz 



Central (C) 


16' 


^59' 


"41f63 


-40° 


03' 


43^'61 


2.0 


± 


1.0 


4.2 


± 


0.8 


8.1 


± 


0.2 


12.1 ±0.2 


Inner-East (1-E) 


16 


59 


41.87 


-40 


03 


44.55 


5.4 


± 


0.9 


7.0 


± 


0.7 


6.3 


± 


0.2 


5.7 ±0.1 


Inner- West (I-W) 


16 


59 


41.35 


-40 


03 


42.94 


5.3 


± 


0.9 


6.7 


± 


0.7 


5.8 


± 


0.2 


3.3 ±0.2 


Outer-East (0-E) 


16 


59 


44.07 


-40 


03 


52.21 


9.7 


± 


0.3 


5.2 


± 


0.2 


9.0 


± 


0.3 


7.6 ±0.4 


Outer- West (0-W) 


16 


59 


39.83 


-40 


03 


41.90 


4.0 


± 


0.7 


8.0 


± 


0.5 


6.2 


± 


0.2 


3.9 ±0.3 



Table 2. SED FITTING 



Parameter Value 



Inclination 


41° 




Stellar Age 


4.5 X 10^ 




Stellar Mass 


14.7 Mq 




Envelope accretion rate 


3.4 X 10" 


3 Mo yr-1 


Disk Mass 


0.26 M0 




Disk accretion rate 


5.5 X 10" 


" Mq yr-i 


Total Luminosity 


5.8 X 10^ 


L0 


Envelope Mass 


1.7 X 10^ 


Mq 
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Fig. 1. — ATCA maps of the radio continuum emission from IRAS 16562-3959. Beams are 
shown in the lower left corner of each panel. Top left: 1.4 GHz map. Contour levels are -4, 
4, 6, 8, 10, 13, 16, 18, and 20 times a {la = 0.35 mJy beam"^). Top right: 2.4 GHz map. 
Contour levels arc -7, 7, 12, 17, 25, 32, 38, 45, 54, and 60 times cj (la = 0.12 mJy beam~^). 
Bottom left: 4.8 GHz map. Contour levels are -4, 3, 6, 9, 13, 19, 28, 40, 60, and 80 times a 
{la — 0.086 mJy beam~^). Bottom right: 8.6 GHz map. Contour levels are -4, 3, 6, 9, 13, 
16, 19, 24, 50, and 100 times a {la = 0.081 mJy beam"^). 
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Fig. 2. — Radio continuum flux density versus frequency for the five principal radio com- 
ponents detected toward IRAS 16562-3959. From top to bottom, spectra of tlie Central, 
Inner-East, Inner- West, Outer-East and Outer- West components. In the top panel the 
dashed hne indicates a least squares power law fit to the data. Empty squares represent 



the data reported by 



Urquhart et al 



(120071 ). In the other panels the dashed line indicates 
a least square fits to the observed spectra with a model of homogeneous region of ionized 
gas. The emission measure (EM) was fitted, and the solid angle {Q) was obtained from the 
deconvolved size of the radio images at 8.6 GHz. Both parameters are given in each of the 

na.np.ls Error bars arp, thp, ma.'X'imiim nf pithp.r thp prrnr va.liips in TablplTlnr 50% 
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Fig. 3. — Gray scale image of the band emission towards IRAS 16562-3959. Black 
contours: 1.2 mm dust continuum emission (SIMBA/SEST). Red contours: 8.6 GHz radio 
continuum emission (ATCA). 



- 23 - 




16'"59'"44= 43' 42' 41' 40' 16''59'"44' 43" 42' 41' 40' 

Right Ascension (.12000) Right Ascension (J2000) 



Fig. 4. — Left Gray image: 8 fim IRAC-Spitzer emission. Black contours: 250 GHz dust 
emission detected with SIMBA. Red contours: 8.6 GHz emission. Crosses mark the position 
of the OH masers. Note that the 8 /im data is saturated near the central source, which 
produces an spurious extension running towards the south. Right: Gray image: band 
emission. Contours and symbols same as in the left panel. 
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Fig. 5. — Grey scale: TIMMI2 10.4 //m emission. Contours: 8.6 GHz radio emission. The 
cross marks the position of the OH maser associated with the central source. There is strong 
10.4 /im emission associated with the central source and more diffuse emission associated 
with the Inner-East lobe. 
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Fig. 6. — Spitzer three color image using the 8.0 fim , 4.5 /im and 3.6 fim IRAC images for 
red, green and blue, respectively. The jet is located at the center of the image associated 
with an extended diffuse "green fuzzy" or enhancement of the 4.5 fim band. 
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Fig. 7. — Spectral energy distribution of the Central source. Crosses: Spitzer/IRAC data 
at 3.6 and 4.5 /im. Empty square: TIMMI2 data at 10.4 fim. Empty circles: MSX data at 
8.3, 12.1, 14.7 and 21.3 jjm. Stars: IRAS data at 12, 25, 60 and 100 jjm . Filled Square: 
SIMBA data at 1200 /xm . The line indicates the best fit using the models described in 



Robitaille et al. 



(120071). 



